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fvq , Abstract. We present results on the evolution of (pt) with energy, for hadrons emerg- 

ing from hadronization of a quark-gluon plasma possibly produced in high-energy 
heavy-ion collisions. We find that BNL-RHIC energy corresponds to the previously 
predicted plateau of (pt), reminiscent of the assumed first-order QCD phase transition. 
Heavy hadrons arc the best messengers of the transverse flow stall. 
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Theoretical arguments [1] and lattice QCD studies [2] suggest that at high tem- 
peratures, roughly T ~ 200 MeV, the thermodynamically stable state of QCD is 
different from that at T = 0, i.e. the theory of strong interactions exhibits a phase 
transition to a new state called "Quark Gluon Plasma" (QGP). It is the only phase 
transition of a fundamental theory accessible to experiments under defined labora- 
tory conditions. In the QGP the effective number of relativistic degrees of freedom 
is expected to be substantially larger than below the phase transition temperature 
T c . The search for signatures of that phase transition is one of the primary goals 
of high-energy heavy-ion collision experiments. 

In particular, if the transition is first order with non- vanishing latent heat, it 
is supposed to leave fingerprints on the hydrodynamical expansion pattern [3,4]. 
Isentropic expansion would proceed through phase coexistence, where the pressure 
p is independent of the energy density e, such that the isentropic speed of sound 
cj. = dp/de = 0. Thus, during that phase coexistence stage energy density gradients 
in the system do not reflect in pressure gradients, and the average transverse flow 
velocity does not increase substantially. As the p t spectra of the emitted hadrons 
are determined by the transverse boost velocity and the temperature at freeze-out, 
one expects a plateau of (p t ) of charged hadrons for some range of beam energies 
where the central region "sweeps" through the mixed phase [3]. 

However, the early studies of that effect suffered from some simplifications which 
render comparisons to existing experimental data, and extrapolation to higher 
energies, unrealistic. First, the latent heat of the transition (and therefore the 
space-time volume of the mixed phase) was overestimated because only pions were 
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assumed to contribute to the entropy of the hadronic phase. Second, final-state 
decays of hadronic resonances (like p — > nn or A — > 7riV) which contribute more 
than 50% to the measured pion multiplicity were not taken into account. Finally, 
the most serious problem arises from the arbitrariness of the decoupling hypersur- 
face, which can not be determined within hydrodynamics but affects the results 
considerably. 

To reduce the uncertainties and the number of free parameters one therefore 
has to improve on the treatment of the post-hadronization stage by considering 
microscopic transport [5,6], e.g. within the relativistic Boltzmann equation with 
collision kernel as determined by the known hadronic cross sections. The expansion 
and subsequent hadronization of the high-temperature QGP state can be modelled 
within relativistic hydrodynamics, assuming that the system evolves locally through 
a phase coexistence. Thus, in the space-time region inbetween the initial-time 
(purely space-like) hypersurface and the hadronization hypersurface (the boundary 
between mixed phase and hadronic phase, which has both space-like and time-like 
parts) we solve the continuity equation for the energy-momentum tensor, 



In the QGP we assume a gas of quasi-free u, d, s quarks and gluons; for sim- 
plicity, the only contribution from interactions to the QGP pressure is due to the 
temperature independent vacuum pressure (bag pressure), see e.g. [1]. 

On the hadronization hypersurface, we switch to the Boltzmann equation, 



where fiix^^p 11 ) denotes the phase-space distribution function of species %. The 
binary collision approximation is used to construct the collision kernel Cj [7], which 
treats almost all hadrons from [8] explicitly. All those hadrons are also taken into 
account in the equation of state of the mixed phase entering eq. (1). The switch 
from hydrodynamics to microscopic transport is performed by matching the energy 
momentum tensors and conserved currents on the hadronization hypersurface. The 
evolution towards freeze-out is thus determined by the relevant hadronic cross- 
sections and decay rates, competing with the local expansion rate d^v? . The initial 
condition for central collisions of Pb/Au nuclei and CERN-SPS energy is that 
hydrodynamic flow sets in on the Tj = 1 fm proper time hypersurface, and the 
entropy per net baryon is s/n B = 45, which reproduces a variety of measured final- 
state hadron multiplicities as well as the p t spectra of ir, K, p, A, S, and fl, see [5,6]. 
For the higher BNL-RHIC energy, ^/s = 200A GeV, the parameters were assumed 
to be Tj = 0.6 fm and s/hb = 205, which yields dN c h/dy ~ 800 at y — 0. Fig. 1 
depicts (p t ) for a variety of hadrons, for CERN-SPS, BNL-RHIC, and CERN-LHC 
energies. One observes that at each energy (p t ) increases with m. Reminiscent 
of the collective expansion before decoupling, all hadrons flow approximately with 
the same transverse velocity, and therefore their momentum essentially increases in 
proportion to their mass. This effect is seen in data obtained for Pb+Pb collisions 
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at the SPS [9]. The spectra of heavy hadrons are much less affected by the random 

thermal motion at decoupling, v t h oc ^jT/m, and are therefore better suited to 
measure the transverse flow build up during the evolution. One also observes that 
despite the five times higher entropy per baryon assumed to be achieved at RH1C 
energy, the average transverse momenta are predicted by this model to be similar 
to the lower SPS energy. This is mainly due to the first-order phase transition 
featuring a region of energy densities where c 2 s is small (it is not exactly zero in this 
calculation due to the finite conserved baryon charge); and to some extent also due 
to earlier decoupling of the hadrons emerging from the hadronization of the plasma 
at RHIC energy [5]. 

If the initial energy density in the central region continues to grow towards higher 
energy as predicted by models assuming saturation of the transverse density of 
gluons [10] (see however [11]), the hydrodynamical expansion will lead to very large 
(pt) at LHC energy, see Fig. 1. Thus, RHIC-energy is possibly right in the plateau 
of (p t ), if the concept of a phase mixture is applicable to high-energy collisions. 

Fig. 2 shows the freeze-out volume of the pions at central rapidity, calculated as 
the average transverse area ir(r t ) 2 times the average length of the central rapidity 
slice, which is equal to proper time r. The space-time volume is essentially in- 
creasing linearly with the multiplicity. Thus, the pion density at freeze-out is not 
affected much by the presence of the phase transition, unlike the average transverse 
momentum. When extrapolating each line to dN^/dy = we obtain a "hollow" 
volume Vo, which increases with p t . This is caused by the radial flow: high p t pions 
can not be emitted from the center, where the collective velocity field vanishes for 
symmetry reasons, and where the pions have only random thermal momenta. 
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FIGURE 1. (pt) versus hadron mass m for n, K, p, A, S, and O. Dots for y/s = 18A GeV, 
squares for y/s = 200A GeV, and triangles for ^/s = 5500A GeV. 
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FIGURE 2. Freeze-out volume of the pions as a function of the pion rapidity density at central 
rapidity, for various pt cuts. Increasingly high p t pions are only emitted from a "shell", the radius 
of the hollow core increasing with pt, as shown in the inset. 



